Nonclassical properties and algebraic characteristics of negative 
binomial states in quantized radiation fields 



Xiao-Guang Wang* Shao-Hua Pan and Guo-Zhen Yang 
CCAST(World Laboratory), P.O. Box 8730, Beijing 100080 
and Laboratory of Optical Physics, Institute of Physics, Chinese Academy of 
Sciences, Beijing 100080, P.R.China 
(February 1, 2008) 

Abstract 

We study the nonclassical properties and algebraic characteristics of the 
negative binomial states introduced by Barnett recently. The ladder opera- 
tor formalism and displacement operator formalism of the negative binomial 
states are found and the algebra involved turns out to be the SU(1,1) Lie 
algebra via the generalized Holstein-Primarkoff realization. These states are 
essentially Peremolov's SU(1,1) coherent states. We reveal their connection 
with the geometric states and find that they are excited geometric states. As 
intermediate states, they interpolate between the number states and geomet- 
ric states. We also point out that they can be recognized as the nonlinear 
coherent states. Their nonclassical properties, such as sub-Poissonian distri- 
bution and squeezing effect are discussed. The quasiprobability distributions 
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in phase space, namely the Q and Wigner functions, are studied in detail. We 

also propose two methods of generation of the negative binomial states. 
PACS numbers:42.50.Dv,03.65.Db,32.80.Pj,42.50.Vk 



2 



Typeset using REVTgX 



I. INTRODUCTION 



Since Stoler et al. introduced the binomial states (BSs)[l], the so-called intermediate 
states have attracted considerable attention of physicists in the field of quantum optics. A 
feature of these states is that they interpolate between two fundamental quantum states, 
such as the number, coherent and squeezed states, and reduce to them in two different 
limits. For instance, the BSs interpolate between the coherent states (the most classi- 
cal) and the number states (the most nonclassical) [1-6] , while the negative binomial states 
(NBSs) interpolate between the coherent states and geometric states[7-ll]. Another feature 
of some intermediate states is that their photon number distributions are some famous dis- 
crete probability distributions in probability theory: the BS corresponds to the binomial 
distribution, the NBS to the negative binomial distribution, the hypergeometric state [12] to 
the hypergeometric distribution, and the negative hypergeometric state[13] to the negative 
hypergeometric distribution. 

Recently Barnett introduced a new definition of NBS [14], 



where \n) is the usual number state, < i] < 1 and M is a no n- negative integer. They 
find that the NBS \q, M) and the BS have similar properties if the roles of the creation 
operator a) and annihilation operator a are interchanged. The photon number probability 
\C n (r), M)\ 2 is associated with the probability that n photons were present given that M 
are found and that the probability for successfully detecting any single photon is r\. Mixed 
states with the photon number probability include those applicable to photodetection and 
optical amplification[15]. 

The BS is a intermediate number- coherent state and the original NBS is a intermediate 
geometric-coherent state. One question naturally arises that if there exist an intermediate 
state which interpolates between the number and geometric state. In reality, the new NBS 



oo 



77, M) = ]T C n (r),M)\n) 



n=M 




(1) 
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is just the intermediate number-geometric state. This fact will be seen in the next sec- 
tion. Thus, we have three intermediate states which interpolate between two of the three 
fundamental states (the number, coherent and geometric states). 

In the present paper we shall study the nonclassical properties and algebraic characteris- 
tics of the new NBS. The ladder operator formalism , displacement operator formalism and 
related algebraic structure will be formulated in section II. It is interesting that the algebraic 
structure is the SU(1,1) Lie algebra via the generalized Holstein-Primarkoff realization. In 
section III, we will show that the NBS can be viewed as excited geometric states, interme- 
diate number-geometric states and nonlinear coherent states. The nonclassical properties, 
such as sub-Poissonian distribution and squeezing effect will be investigated in detail in 
section IV. The Q and Wigner functions are studied in Section V and the two methods of 
generation of the NBS are proposed in section VI. A conclusion is given in section VII. 

II. LADDER OPERATOR FORMALISM DISPLACEMENT OPERATOR 
FORMALISM AND ALGEBRAIC STRUCTURE OF THE NEW NBS 

A. Ladder operator formalism and algebraic structure 

It is known that the BSs are special SU(2) coherent states[5,6] and the original NBSs are 
Perelomov's SU(1,1) coherent states[ll] via the standard Holstein-Primakoff realizations. 
So we expect that the algebra involved in the new NBS is SU(1,1) Lie algebra. 

It is easy to evaluate that 



a tn \r],M) 



(M + n)\ 



1/2 



\r],M + n). (2) 



M\r] n 

In particular, for n = 1, we get 
/ M + 1 \ 1/2 

a%M) = (— ^) \V,M + 1). (3) 

The creation operator raises the NBS \rj : M) to \rj : M + 1). This property is similar to the 
action of the creation operator on the Fock state |M), 



(J\M) = v / MTT|M + l). (4) 



Actually, in the limit of rj — > 1, the NBS \rj,M) reduces to the number state \M) and 
Eq.(3) naturally reduces to Eq.(4). 



The key and interesting point is that there exists another operator y N — M which also 
raises the NBS \rj : M) to \r/, M + 1). From Eq.(l), the following equation is directly derived 



as 



^n-m\t],m) = {^-^j VmTT\v,m + 1). (5) 

Comparing Eq.(3) and Eq.(5), we get 



\In-M\ti, M) = yjl - rja^rj, M). (6) 

Multiplying the both sides of the above equation by the operator \J N — M from left, we 
obtain the ladder operator formalism of the NBS as 

(N-Jl- r]\/N - M<f)\ri, M) = M\rj, M). (7) 

In the limit of 77 — > 1, we find that the NBS|r], M) reduces to the number state \M) and 
Eq.(7) to the equation N\M) = M\M) as expected. 

It can be proved that the operators appearing in Eq.(7) can form SU(1,1) Lie algebra: 

K = N - M ~ - , K + = \f N — Ma\ K_ = a\J N — M. (8) 
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Reminding of the standard Holstein-Primarkoff realization of SU(1,1) Lie algebra, J< 







N + ^-, J + = y M + N — la\ J_ = ay M + N — 1, we call the new realization as generalized 
Holstein-Primakoff realization. To our knowledge, the new realization of SU(1,1) Lie algebra 
seems not to be addressed in the literature. 

In terms of the generators K Q , K + and K_ of the SU(1,1) Lie algebra , Eq.(7) is rewritten 

as 

(K - y^K + )\ V ,M) = M±l\r,,M). (9) 
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B. Displacement operator formalism 

Now we try to find the displacement operator formalism of the NBS. To this end, let us 
rewrite the NBS as 

~ 'M + n} 



oo 

M+l 



n=0 

oo 



1/2 



\M + n) 



= ^E (Vr 74 atr iM>. do) 



n=0 VH! 

Then by making use of the following identity 

[f(N)a*] n = a^ n f(N + n)f(N + n- l)...f(N + 1), (11) 
Eq.(10) can be written in the exponential form 

\ V ,M) = r ] M ^e^ K +\M}. (12) 

Note that K-\M) = 0, then the displacement operator formalism is obtained as 

\r),M) = e VT= ^ K+ 'q Ko e-' /T ^ K -\M) 

= e^ K+ - K -^\M), (13) 

where £=arctanh y/1 — rj. In the derivation of the above equation, we have used the identity 

e aK + - a -K- = e7 ^ +(1 _ | 7 | 2) ^o e -7^- 5 ( 14 ) 

where 7 = atanh j ce j / j a; j . As seen from Eq.(13), the NBS can be simply recognized as 
SU(1,1) displaced number states. Actually, the NBS are essentially Peremolov's coherent 
states as shown below. 
On the space 

S = span{|n + M) = \n; k)\n — 0, 1, 2...} , k — (15) 
we have 



K + \n;k) = \J(n + l)(2k + n)\n + 1; k) 



K_ \n; k) = ^Jn(2k + n- l)\n - 1; k) 

K \n;k) = (n + k)\n;k) (16) 



This is the discrete representation of SU(1,1) Lie algebra with Bargaman index k = 
(M + l)/2. We see that the generalized Holstein-Primakoff realization gives rise to the 
representation of SU(1,1) on the space S. Note that \M) = |0; k), the NBS can be written 

as 

\ V ,M) = e ^ K +- K - ) \0;k). (17) 
This shows the NBS are essentially Peremolov's coherent states. 

III. THE NBS AS EXCITED GEOMETRIC STATE, INTERMEDIATE 
NUMBER-GEOMETRIC STATE AND NONLINEAR COHERENT STATE 

A. As excited geometric states and intermediate number-geometric states 

From Eq.(lO) we obtain 



'M + n x 1/2 



n=0 



n) 



(M+l)/2 oo 

V M ! n=0 

M/2 

" r« tM !^ ( 18 ) 



VM! 
where 

oo 

\V) 9 = V 1,2 T.^-V) n ' 2 \n) (19) 

ra=0 

is the geometric state [16-21], which is also called Susskind-Glogower phase state [11], phase 
eigenstate[16,17],and coherent phase state[19]. The photon number distribution is 77(1 — i]) n , 
the geometric distribution. From Eq.(18) the NBSs can be generated by repeated application 
of the creation operator at on the geometric states. This shows that the NBS belongs to 
an interesting class of nonclassical states, excited quantum states. These states are first 
introduced by Agarwal and Tara as excited coherent states[22]. So the NBSs can be viewed 
as excited geometric states. 
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Setting M = in Eq.(18), we get 

faO) = \v) a , (20) 

which shows that the NBS reduces to the geometric states for M — 0. Note that the NBS 
reduces to the number state \M) in the limit of i] — > 1. Thus, the NBS interpolates between 
the number state and geometric state and can be viewed as number-geometric state. 

B. As nonlinear coherent states 

The geometric states are the eigenstates of the Susskind-Glogower phase operator[23] 
(1 + N) l l 2 a, obeying the equation 

{N + \)-V 2 a\n) a = y[r^i\r,) r (21) 

Comparing with the definition of the nonlinear coherent states|a) n ; [24,25] 

f(N)a\a) nl = a\a) nh (22) 

we know that the geometric states are nonlinear coherent states with the nonlinear function 
f(N) = (N + l)-W 

In a previous work, we have proved a general result that the excited nonlinear coherent 
states are still nonlinear coherent states[26]. Since the geometric states are nonlinear coher- 
ent states and the NBSs can be recognized as excited geometric states, we infer that the 
NBSs are nonlinear coherent states. In fact, multiplying Eq.(7) by the annihilation operator 
a from the left , we get 



(N + 1 - M)a\q, M) = yjl- r]VN + 1 - M(l + N)\q, M). (23) 

Since we discuss the problem in the space S^Eq^lS)), we can multiply Eq.(23) by 1/[(1 + 
N)y/N + 1 - M] from left. This leads to 

[VN + 1 - M/(N + l)]a|?7, M) = yjl- rj\v, M). (24) 

The above equation shows that the NBS are nonlinear coherent states with the nonlinear 
function + 1 - M/(N + 1). Eq.(24) naturally reduces to Eq.(21) for M = 0. 



IV. NONCLASSICAL PROPERTIES 



A. Sub-Poissonian distribution 

The simplest way to investigate the statistical characteristics of the radiation field is to 
differentiate the generation function 

G( A)= S i.( B )V = ^( TT A- I ) (25) 

with respect to the auxiliary real number A. Here P{n) = \C n (r), M)\ 2 is the photon distri- 
bution function of the NBS. The factorial moments are defined as F(n) = ^^|a=i- From 
Eq.(25) we obtain the factorial moments F(l) and F(2) as 

F(l) = (N) = ^±1 - l, (26) 

F(2) = (J V> - (N) = (M + 2 >< M+1) - 4^±I + 2. (27) 

I] 2 T] 

Then we can easily derive Mandel's Q parameter 
(N 2 ) - (N) 2 - (N) 



Q 



(N) 

F(2)-F 2 (l) r] 2 -2(M +l)r] + M +1 



(28) 



F(l) 77(M+1-t7) 
which measures the deviation from the Poisson distribution which corresponds to the co- 
herent state with Q — 0. If Q < 0(> 0), the field is called sub(super)-Poissonian. The 
denominator of Eq.(28) is positive since rj < 1, while the numerator can be positive or neg- 
ative. For M — 0, Q — (1 — 77) / 77 > 0. The NBS (77, 0) (geometric state) is super-Poissonian 
except 77 = 1. For M > 0, the condition for the numerator r? 2 — 2r/(M + l) + M + l<0is 

v > = M + 1 - ^M(M + 1). (29) 

It can be proved that < ?7_ < 1. Thus the NBS \rj,M)~(M > 0) is super-Poissonian 
when 77 < t7_ and sub-Poissonian when 77 > 77 . As M increases, 77_ decreases and the 
sub-Poissonian range increases. 
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B. Squeezing effect 



Define the quadrature operators X (coordinate) and Y (momentum) by 

X=^(a + at),F = i(a-at). (30) 
Then their variances 

Var(X) = (X 2 ) - {X} 2 , (Var(Y) = (Y 2 ) - (Y) 2 (31) 

obey the Heisenberg's uncertainty relation 

Var(X)Var(Y) > (32) 

If one of the Var(X) and Var(Y) is less than 1/4, the squeezing occurs. In the present case, 
(a) and (a 2 ) are real. Thus, the variances of X and Y can be written as 

Var(X) = I + l« a ta) + (a 2 ) - 2(a) 2 ), (33) 
Var(Y) = \ + \{{<Ja)-{a 2 )). (34) 

From Eq.(2), the expectation values (a) and (a 2 ) are obtained as 

(a) = £ f^^'Y^V'V-^^ + l) 172 (35) 



n=M 



(a 2 ) = r, M +\l- V )- M ± ( n + 2 y /2 ^ 1/2 (l-^r 1 [(- + 2)(n+l)] 1 / 2 (36) 

Using Eqs.(26), and (33)-(36), we can investigate the squeezing effect. 

By numerical calculations, we find that the squeezing occurs in both the quadrature X 
and Y. Fig.l gives the variance of the quadrature X versus rj for different M. It can be seen 
that the range and degree of squeezing increase as M increases. There exists a critical value 
of M. When M < 7, there is no squeezing for arbitrary values of rj. The squeezing also 
occurs in the quadrature Y as shown in Fig. 2. In contrary to the squeezing in the quadrature 
X, the range and degree decrease as M increases. When M is larger than a critical value 
31, no squeezing occurs. 
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V. THE Q AND WIGNER FUNCTIONS 



If a field is prepared in a quantum state described by a density operator p, we can define 
the s-parametrized quasiprobability distribution in phase space as[27] 

P((3,s) = ^J d 2 iC{i- s) exp(/3f - (3*0 , (37) 

where the quantum characteristic function is 

C(Z,s) = Tr[D(Op}exp(s\Z\ 2 /2). (38) 

Here (5 = x + iy, with (x, y) being the c numbers corresponding to the quadratures (X, Y), 
and -D(£) = exp^a 1 " — £*a) is Glauber's displacement operator. It is possible to write the 
s-parametrized quasiprobability distribution as a infinite series [28] 
2 00 ( 1 4- <i) k 

P(P,s) = -E(-i) n_ (39) 

71 n=0 S ^ 

where fc) = -D(/3)]A;) is the so-called displaced number state. The expression above is 
suitable for direct numerical calculations. 



A. Q function 

If we take s = — 1, Eq.(39) is reduced to the familiar expression for the Q function 

Q{(3) = -((3\p\(3). (40) 

71 

Note that since the NBS can be viewed as an excited geometric state [see Eq.(18)], we obtain 
the Q function of the NBS as 

2 

/Ml. (41) 

In Fig. 3 we present plots of the Q function of a NBS for different values of rj and M = 5. We 
can clearly see the deformation of the Q function. When rj — 1, the Q function representing 
the number state \M — 5) is formed[Fig3.(d)] as expected. From the Q function we can 
also study the squeezing effects by examining the deformation of their contours. Fig. 4 is the 
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Q((3)= V M+1 exp(- 



2M 



J2(3 n (l-r]) n / 2 /Vri. 



n=0 



contour plot of the Q functions of two particular NBSs, with (a) r\ = 0.2, M = 5 and (b) 
rj = 0.85, M = 50. We can clearly see the compression along the y and x direction, which 
corresponds to squeezing in the Y quadrature of the first NBS and in the X quadrature of 
the second. 



B. Wigner function 

By taking s = in Eq.(39), we obtain a series representation for the Wigner function 

r\ OO 

W(P) = -J2(-V k (P,k\p\[3,k). (42) 

, n 

k=0 



Now we insert Eq.(l) into the above equation, which yields 



W{(3) = -j: 



2 <x x ~ 



71 k=0 



J2 C n ( V ,M)Xnk(P) 



n=M 



(43) 



In the expression above, the matrix elements Xnk(P) = (n\D((3)\k) are given by[29] 



Xnk(P) = P n (-Pl k eM-\P\V^F (-n,-k;\P\- 2 )/V^., (44) 

where 2 Fo(a, P; z) are the generalized hypergeometric functions [30]. The present form of 
Xnk(P) is convenient for numerical calculations. 

The Wigner function can be used to trace the nonclassical behaviors of quantum states. 
It is known that the negativity of the Wigner function is a sufficient but not necessary 
condition for having nonclassical effects. In Fig. 5 we give plots of the Wigner function of 
a NBS by numerical calculations of Eq.(43) for different values of i] and M — 1. As in 
Fig. 5(a), the negative part of the Wigner function is already noticeable for rj = 0.3. For 
rj = 0.5[Fig.5(b)], the negative part is pronounced. In fig5(c), for rj = 0.9, the negative part 
is even larger, and finally, in Fig. 7(d), we have the full Wigner function of a number state 
|1)(?7 = 1). The Wigner function becomes more and more negative as 77 increases. 
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VI. GENERATION OF THE NEW NBS 



Let us discuss the dynamical generation of the NBS. We consider two different meth- 
ods. The first is quite straightforward in concept from the displacement operator formal- 
ism(Eq.(13)) but might not be very easy to achieve experimentally. The Hamiltonian is 
given by 



H = H + ixl^N - Mat exp(-iwt) - a\J N - M exp(icut)] 

H = ua)a. (45) 

The constant x is t ne coupling strength. The coupling is of intensity-dependent type and is 
similar to those in some intensity-dependent Jaynes-Cummings models [31-34]. The unitary 
time evolution operator in the interaction picture is 

U(t) = exp[xt(V 'N - Ma f - cl\Jn - M)]. (46) 

Supposing the system is initially prepared in the number state \M), we find the system 
at time t is the NBS 

U(t) | M) = 1 1 - tanh 2 ( X t) , M) . (47) 

The second method of the generation of the NBS is based on the fact that the NBS is the 
excited geometric state. The geometric state can be prepared in the non-degenerate three- 
wave interaction system[35]. We can also generate the geometric state by the non-degenerate 
parametric amplifier described by the two-mode Hamiltonian [36] 

H = H + ix[a\a\ exp(— 2iut) — aia 2 exp(2ia;£)], 

H = uj x a\ai + uj 2 a\a 2 , (48) 

where a\ and cu 1 (a 2 and uj 2 ) are the annihilation operator and frequency for the signal(idler) 
mode. Frequencies ui\ and u 2 sum to the pump frequency, 2oj = oj\ + uj 2 . The coupling 
constant x is proportional to the second-order susceptibility of the medium and to the 
amplitude of the pump. The unitary time evolution operator in the interaction picture is 

U(t) = exp[xt(a\al- ai a 2 )] (49) 
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Suppose that the system is initially prepared in the state |0,0) = |0)i <E> |0) 2 . Then at 
any time t the system is in the state 

oo 

\r ] } tm = U(t)\0,0}=r ] 1 / 2 J2(±-vT /2 \n,n) (50) 

n=0 

which is the two-mode geometric state in comparison with Eq.(19). Here rj = 1 — tanh 2 (%t). 

Once the two-mode geometric state is prepared, one can generate two-mode NBS by 
the following procedure in analogy to that proposed by Agawal and Tara[22]. Consider the 
passage of a two-level excited atom through a cavity. Let the initial state of the atom-field 
system be \r}) tm <8> |e), where |e) is the atomic excited state. The interaction Hamiltonian 
has the form [3 7] 

H = h(gS + ai + g*S~a{), (51) 

where S' ± are the psedospin operators of the atom and g is the coupling constant. Since g 
is generally small, the state at time t can be approximated by 

|V>(0> ~ \v)tm ® \e) - ig*ta\\rj) tm ® \g), (52) 

which is valid for interaction times gt << 1. From the above equation we observe that, if the 
atom is detected to be in the ground state \g), then the state of the field is reduced to a\ \rj) tm . 
An extension of the above arguments to the multiphoton processes would imply that the 
state a\ M \r])t m can be produced in multiphoton processes. For the multiphoton processes, 
the Hamiltonian (Eq.(51)) is replaced by a new Hamiltonian with a x af 1 . Thus, the above 
procedure for a multiphoton process will result in the state 

fa M) tm = V {M+m £ H ) V V " VT /2 \M + n,n). (53) 

The above normalized state is just the two- mode NBS. In a short summary two methods 
are proposed to generate the NBS. 
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VII. CONCLUSIONS 



We have investigated the NBS induced recently by Barnett and found the ladder operator 
formalism and displacement operator formalism of the NBS. The algebra involved is SU(1,1) 
Lie algebra via the generalized Holstein-Primakoff realization. We found that the NBS are 
essentially Peremolov's coherent states. 

As excited quantum states, the NBSs are excited geometric states. As intermediate 
states, they interpolate between the number and geometric states. According to the defini- 
tion of the nonlinear coherent states, we find that the NBSs are nonlinear coherent states. 

The NBS can be sub-Poissonian or super-Poissonian. There exists a critical point rj-. 
The NBS is sub-Poissonian when rj > rj- and super-Poissonian when rj < rj-. The squeezing 
occurs in both the quadrature X and Y with two critical values of M = 7 and 31, respectively. 
There is no squeezing occurs in the quadrature X for M < 7 and in the quadrature Y for 
M > 31 for arbitrary values of rj. The Q and Wigner functions are studied numerically. 
They show that the NBSs have prominent nonclassical properties. We have proposed two 
methods of generation of the new NBS. 

In addition, the remarkable properties of the new NBS seem to suggest that it deserves 
further attention from both theoretical and application sides of quantum optics. 
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FIG. 1. The variance of the quatrature X as a function r\ for different values of M: (a)M = 5, (b)M = 10, and 
(c)M = 35. 
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Fig. 2. The variance of the quatrature Y as a function r\ for different values of M: (a)M = 5, (b)M = 20, 
and (c)M = 40. 
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Fig. 3. Q function of a NBS for four values of rj: (a)rj = 0.3, (b)rj = 0.5, (c)rj = 0.9, and (d)rj = 1. In all 
cases M = 5. 
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Fig.4 Contours of the Q function of two NBSs with (a)ry = 0.2, M = 5, and (b)»j = 0.85, M = 50 
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Fig. 5. Wigner function of a NBS for four values of rj: (a)rj = 0.3, (h)rj = 0.5, (c)rj = 0.9, and (d)r] = 1. In 
all cases M = 1. 
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